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The dynamic response functions of strongly interacting fermion gas in homogeneous space are 
investigated in a virial expansion to second order. The density response function exhibits transition 
from atomic to molecular response, as the interaction strength increases and the system undergoes 
BCS-BEC crossover. The qualitative features of density and spin response agree with measurements 
from the Bragg spectroscopy experiments. The virial response is exact at low density and high 
temperature, therefore providing a benchmark for many-body response. 
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I. INTRODUCTION 

Recent years the advances in cold atom experiments have improved our understanding of thermodynamic and 
dynamic properties of strongly interacting systems. Particularly, the dynamic response function gives the response 
of the system under an excitation probe, therefore describing the inelastic scattering processes. Quite remarkably, 
the Bragg spectroscopy has been used to measure the dynamic and static density response functions for a strongly 
interacting cold Fermion system ST. A smooth transition from atomic to molecular spectra, or from BCS to BEC 
regimes, is observed with a clear signature of pairing at and above unitarity, which provides a direct link to two-body 
correlations. Recently, the Bragg spectroscopy is also applied to obtain the dynamic spin response function [2J. 

Due to its nonperturbative nature, theoretically it is extremely hard to obtain the response function of strongly 
interacting Fermion system. There have been several investigations on the dynamic response function [3]-[7] . In certain 
limits, large momentum transfer between external probe and the system [8], or system being at very low density, to 
name a few, the many-body phenomena are expected to be dominated by few-body physics explicitly. The virial 
expansion presents a tractable approach to the strongly interacting system and has a controllable small parameter, 
the fugacity z = exp(/i/T), when the system is at low density (chemical potential /i) and high temperature (T). It has 
been applied to study the thermodynamic properties of strongly interacting Fermi gas (9j 110) . For a trapped strongly 
interacting Fermi gas, a quantum virial expansion to second order for the response functions has been developed in 

Refs. puna. 

In this work, the latter virial expansion is extended to study the dynamic response functions of strongly interacting 
Fermi gas in homogeneous space. The formalism closely follows that developed in Ref. [TT] for trapped Fermi gas. 
The dynamic density response function in homogeneous space is found to exhibit transition from atomic to molecular 
spectra, as the interaction strength increases and the system undergoes BCS-BEC crossover. The virial response is 
exact at low density and high temperature, when fugacity is a small number. Therefore it provides a benchmark for the 
many-body response functions. Qualitatively the response functions of strongly interacting fermi gas in homegeneous 
space show similar characteristics as those in trapped fermi gas. In fact they may be related by a local density 
approximation. In this work, the results in homogeneous space are explicitly written down in compact and closed 
form with simple integrals. The virial expansion for dynamic response function can be readily applied to other strongly 
interacting many-body system, like neutron matter and nuclear matter in supernova (see Ref. [13 for a study on the 
long wavelength behavior of static response function). 

The paper is organized as follows. In next section [TT] the formalism for the dynamic response function of Fermi gas 
in second order virial expansion is presented in detail. In Sec. [ill] the frequency dependence, interaction dependence, 
and temperature dependence of dynamic response function, as well as the static response function are discussed. The 
conclusions are given in Sec. |IV| 

The natural units h = ks = 1 are used throughout. 
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II. THE DYNAMIC DENSITY RESPONSE 

The formalism below closely follows that developed in Ref. [11] for a trapped Fermi gas. For completeness the 
formalism is given in detail here. The dynamic density response for a spin-unpolarized gas under four momentum 
transfer (q,w) can be separated into two pieces: 

S D (< l ,u)=2[S n (f l ,u>) + S n (q,u)], (2.1) 

where = and = are used. Similarly one could obtain the dynamic spin response function, 

S s (q,w) = 2[S rt (q,w) - S n (q,w)]. (2.2) 

It is most convenient to relate the response function to a time-dependent correlation function in imaginary time, via 
analytic continuation and the fluctuation-dissipation theorem, i.e. the dynamic susceptibility, 

W(r,r» = -<T T 7v(r,T)rv(r,0)) (2.3) 

where T T is the imaginary-time ordering operator, and ?v(r, r) is the density (fluctuation) operator in spin channel 
cr, and r is an imaginary time in the interval < r < (3 = 1/ksT. With these definitions, one can show that 

q (r r > : A - ~ Im ^' ( r ' r ' ; itJn ^ w + * 0+ ) fo A\ 

*aa'( r » r = W )- - e _^j . ( 2 - 4 ) 

where the matsubara frequencies w„ = 2imkBT (n — 0,±1,...). Above equation can be Fourier-transformed with 
respect to r - r' to obtain S acr ,(q,u). By the usual virial expansion in fugacity z — e^ M , 

W = zX 1 [l + z(X 2 /X 1 -Q 1 ) + ...}, (2.5) 

where 



X N — -Tr N 



e-^e^n CT (r)e-^n^(r') . (2.6) 



The interaction enters from second order term. Using completeness relations, 

Ax CTCT ,, 2 (r,r',r) = AX 2 =X 2 - A 2 ° = (2.7) 

= -E [e-^ +T{E ^ ] C P J{v,v')\ , (2.8) 
P,Q 

where 

C P J{v,v>) = {P\n a {v)\Q){Q\h a '{v')\P) (2.9) 

and pair state (tl) \P) (\Q)) has energy E P (Eg). One can perform discrete Fourier transform on Ax^a'pi*, r',r), 
and analytical continuate the results via iw„ — > u; + i?/, to get the second order response function, 

A^w, 2 (r,r, W ) = (1 _ e _^ )7r 

P,Q 

The remaining step to get momentum space response function is through Fourier transform as follows, 

AS CT(T%2 (q,w) = J dRe- iQ ° R J drdT'e ir ^ 2 -^ +ir '^ Q "/ 2+ ^AS^ t2 (T,T',uj) (2.11) 

= ^^ + J Bp- J B Q )e-^F; CT < ?(q), (2.12) 
PQ 
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where 

iC?(q) = J dvdv'e-^-^C P J,{v,v'). (2.13) 
A reminder here for the form of C^®(r, r'), 

< T Q (r,r') = ydr 2 dr' 2 $^(r,r 2 )$ Q (r,r 2 )$^(r',^)$ P (r',r^), 

C t P 4 Q (r,r') - J 'dr 1 dr 2 $^(r,r a )$o(r,r a )*^(r 1 ,r , )$p(r 1) r') ) (2.14) 
One can further decompose the response function in center of mass (CM.) and relative parts, by making use of 

5{lo + E P - E Q )e- f3Ep = j dJ5{J + € pl -t ql )8{u-J + e p2 -eq2)e- !3{t ^ + ^\ (2.15) 

where e p i( g i), e p 2( 9 2) are the CM. and relative part of energy. 

Now one has to differentiate between and t-i cases. For f1\ one needs to group (r, r 2 ), and (r',r' 2 ) separately. 
Spcfically, define CM and relative parts as follows, 

r = R + x/2, 

r 2 = R-x/2, (2.16) 

and do the same for (r'r 2 ). For t-ii one needs to group (r,r 2 ), and (i"i,r') separately. After some algebra, one can 
arrive at following forms, 

= |y dw , W cm (q,w , )TO'(q,w- w')} , (2-17) 



with following notations, 



where, 



W cm = ]T + epi - e 9 i)e"^ pl | 2 , (2.18) 



/ pl9l = J dKe-^^iR^R), (2.19) 



and 



e^^- W ' + e p2 - e^e"^^, (2.20) 

p2q2 

with A^ g2 = \A p2q2 \ 2 and A^ ?2 = A p2 , 2 ,4 ?2 p 2 , and 

A p2q2 = J dxe-^- x / 2 ( /»; 2 (x) ( /» 92 (x). (2.21) 

The center of mass piece is trivial, 

W cm {J) = J^e-^'-e/trnfPrn/cf^ (2 22) 

nir z pq 

where the number of particles has been divided. Next the relative piece is discussed in detail, from BCS to BEC side. 
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A. BCS side 



Using the scattering state wave function <j>sipx) in Eq. (Al |, the relative piece in Eq. (2.201 can be written as, 

W? e i(q,0j) = 



dp 2 p\dq 2 q\8[uj + e P2 - e q2 )e 0€p 



dxx jo(qx/2)4>s(P2x)(j)s(q2x) 

\ 2 



dxx ji(qx/2)\j -ji(p 2 x)(j)s(q2x) 



l>0 \ J 

+ I J dxx 2 ji{qx/2)^ji{q 2 x)<j)s{p2x)\ - non. inter, terms 



(2.23) 



Note either 4>s{p2x) or cf>s(q2x) has to be interacting 5 wave state, otherwise it will cancel with corresponding non- 
interacting term. With above equation and Eq. (2.22) one can proceed to the dynamic structure function in Eq. 



# 2 p^g 2 ^^ e -("+2^-2gJ- 9 -V2) 2 /2? 2 f-2p=/f 



dxx jo(qx/2)cf) S (p 2 x)(j)s(q 2 x) 

2 



^(2/ + l)(-l)'( 1 -^') ( f dxx 2 n {qx/2)J- n {p2x)4>s(q2x) 
i>o V 

^2(21 + 1)(— ' } ( J dxx 2 ji{qx/2)^ji(q 2 x)(t>s(p2x)\ - non. inter, terms 



(2.24) 



where scaled variables, momentum in terms of kp and energy in terms of ep = k F /2m, have been used. For instance, 
q = q/kp, f = T/e F , and AS aa ' j2 = A5 ffff - i2 x e F . 



a) The first term, S wave to S wave scattering, (and its non-interacting counterpart) in Eq. (2.24) can be worked 
out as follows, 



d p 2d q 2 ^(lf e -(^pl-^'fm 2 nff-2pl/f x 

-C p2 C q2 (-sign(q/2 - p 2 - q 2 ) + sign(g/2 + p 2 - q 2 ) + sign(q/2-p 2 + q 2 ) - sign(g/2 + p 2 + q 2 )) 
+ ^S P 2Sq2{s\gn{q/2 -p2- q 2 ) + sign(q/2 + p 2 - q 2 ) + sign(<?/2 - p 2 + q 2 ) + sign(q/2+p 2 + q 2 )) 



C p2 S q2 log 



(l-v 2 f 



C q2 S p2 log 



(l-« 3 ) 2 



— non. inter, terms, 



(2.25) 



where sign() is the sign function, and following quantities are defined, 



C p 2 — 


1 


(2.26) 




Sp2 — 


-p 2 kFa 


(2.27) 


y/l+p^a 2 ' 


V2 = 


q 2 /4 + p 2 -q 2 


(2.28) 


qp2 


V3 = 


q 2 /4 + ql-p 2 2 
q~q~2 


(2.29) 



b) The second and third terms in Eq. ( 2.24 ) due to transition between I > partial waves and S wave are simplified 
as follows, 



dp 2 dq 2 -—r{-) e 
Aq 6 7r 



J, ■ici+2p|-2q=-g 2 /2) 2 /2? 2 T-2p|/f ^(2Z + 1) (- 1)^ 1 "^' ) ; 

1>Q 
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(g ) 2p iM 2 (C 2 q2 - 1) + QtfafSfa + 7TP i (V2)Qi(^)C 8 2^2 

+ (lfPiM 2 (C 2 p2 - 1) + Q ; (« 3 ) 2 ^ 2 + 7rP ; (^)Q^3)C p2 S p2 



(2.30) 



where p and Qi are Z*' l -order Legendre polynomials of first and second kind, respectively. Note except (Qi) 2 terms, 
the other terms are nonzero only for Iuq] or \vs\ < 1 |14) . Using the summation relation for the product of two 
Legendre polynomials (see Appendix |Bj, b) term can be further simplified as follows. 

1. Spin anti-parallel case 



Aq A 7r 

C-) 2 {5{v 2 ) - l)(C 2 q2 1) + {±- log(^) 2 - Qo(v 2 ) 2 )S 2 q2 + Q (v 2 ))C q2 S q2 

+ (l) 2 (S(v 3 ) - l)(C 2 p2 - 1) + (-L log(^) 2 - Q (v 3 ) 2 )S 2 2 + tt(~ - Q (v 3 ))C p2 S p2 . 

(2.31) 

Note again the 1st, 3rd, 4th, and 6th terms are nonzero only for \v 2 \ or |«3| < 1. Principal-vaule integrals are 
assumed where it is appropriate. 

2. Spin parallel case 



AS| T2 (q,tj) 



4g J 7T 



-(^) 2 PoM(C 2 q2 - 1) + (-5— - Qo(v2) 2 )S 2 q2 - irPo(v 2 )Q Q (v 2 )C q2 S q2 
i v 2 i 

<l) 2 PoM(C 2 p2 - 1) + - Qo(f 3 ) 2 )S 2 2 - 7rPo(v 3 )Qo(«3)C p2 5 p2 



(2.32) 



Similarly, the 1st, 3rd, 4th, and 6th terms are nonzero only for \v 2 \ or \v 3 \ < 1. 



Summing up Eqs. (2.25) with (2.311 or (2.32 1, one could obtain the dynamic response functions on the BCS side as 
follows, 



Sff s (q,u) = Z 2 (A^, 2 (g, W )+A^ 2 ( g , W )), 
5 t s T cs (g, W ) = § P (q,u) + z 2 (AS^ 2 (q,u) + AS^q^)), 

where Sp(q,u>) is the dynamic response function for non-interacting Fermi gas [15j . 



(2.33) 
(2.34) 



S F (q,u) 



3T 



1 



IQq l - e-"/ T 



log 



1 + ze 



-(u/q-q) 2 /iT 



1 + ze-(' i /9+'?) 2 /4r 



Note a factor of half is included in above equation to remove spin degeneracy. 



B. BEC side 



(2.35) 



On BEC side there exist contributions from bound state (j>b{x) in Eq. (A2|, besides scattering states, 
c) Transition between bound state and scattering state 



d q 2 q2^ e -(u>-2/k%a 2 -2fi-f/2) 2 /2q 2 f+2/fk 2 F a 

2 4g 

£(2Z + i)(-r/d- 

;>o 

/ dp 2 p 2 — e 



dxx 2 j (qx/2)(j) S (q 2 x) VAtt^ (x) 



dxx ji(qx/2)\ -ji(q 2 x)vAir4>b(x) 



.-. ~2 „-{Cb+2p%+2/k 2 F a 2 -q 2 /2) 2 /2q 2 f-2pl/f 



dxx 2 ji(qx/2)cf) S (p 2 x)V 4TT(j) b (x) 



G 



^Z + lX-l)^- 5 -') \J dxx 2 j l {qx/2) S j'^j l {p 2 x)^()> b {x) 

~ 3 T -(Q-2/k%a 2 -2qZ-q 2 /2) 2 /2q 2 f+2/fkla 2 

?2 4«f 

+ {Qo{u 2 ) 2 C q2 + (arctan[(g/2 - q2)k F a] + arctan[(g/2 + ^kpa^S^ 
3T 



{A(u 2 ) - Q {u 2 f 



(A(u 3 ) - Q0M 2 



+ (Qo(u 3 ) 2 C P2 + (&rcta,n[(q/2 - p 2 )k F a\ + arctan[(q/2 + p 2 )kpa])S P2 ] 
where the following notations have been used, 

l/k 2 F a 2 + q 2 /A + q 2 



U2 = 



U3 



l/k 2 F a 2 + g 2 /4 + p| 



qp2 



A(x) 
A{x) 



1 



x 2 - V 
Qo(l/x)/x, 



for spin ft, 
for spin f| ■ 



(2.36) 

(2.37) 
(2.38) 

(2.39) 
(2.40) 



d) Transition between bound states, or molecular response 



AS*, 2 {q,u) = ?L e -(*~?/lf/2?f + 2/f k 2 F a 2 
Aq 



dxx 2 jo {qx/2)A-K(j)b{xY 



3ttT 
2q 2 kpa 



- { *-q 2 /2) 2 ,2 q 2 f + 2,fkW {l + - 2fc 2 2)-l/2p-V2 [(1 + ^2^-1/ 



(2.41) 



It is clear from Eq. (2.41) that on deep BEC side, kpa <C 1, the molecular response dominates in the dynamic 
response function. 



functions on the BEC side as follows, 



Summing up Eqs. (2.25) with (2.31) or (2.32), as well as (2.36) and (2.41), one could obtain the dynamic response 



Sff c (q,u) = z\AS^ 2 (q,cj) + A% 3 ( 5)W ) + AS^ 2 (q, w) + AS^ 2 (q,u>)), 



(2.42) 



S^ EC (q,uj) = S F (q,u) + z 2 (AS^ t2 (q,ui)+AS^ t2 (q,Lj) + AS^ :2 {q,u;)+AS^ t2 (q,u;)), (2.43) 

III. RESULTS AND DISCUSSIONS 

In this section, the dependences of dynamic response function on frequency, interaction, and temperature, as well 
as static response function are discussed in details. The Bragg spectroscopy experiments have been carried out using 
large momentum transfer. Therefore in this section, the response functions are calculated with a large momentum 
transfer q = 3k p. 



A. Frequency dependence of Sff and Sf± 

In Figure [Tl the dynamic response functions for spin-parallel and spin-anti-parallel cases, and S 1 -^, are shown 
as function of frequency at different interaction strengths for T = 2Tp (left panels) and T = 3Tp (right panels). In 
all the results, the largest fugazity is z—0.25 (BCS side \jkpa = —1 and T/Tp—2). These small values suggest that 
higher order corrections to the virial expansion will be small, ton = q 2 /2m is the atomic recoil frequency where m is 
the atomic mass. 

The effect of interaction is most clear in the dynamic response function for spin-anti-parallel case, S^, as Eqs. 
(2.33 2.42) contain only interaction induced contributions. On the BCS side 1/kpa = — 1, S^i is peaked around the 
molecular recoil frequency u>r/2 due to strong pair correlation from attractive interaction. The peak becomes more 
visible as temperature decreases and/or interactions strength increases. On the BEC side, the response is sharply 
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function 
Eq, 



peaked around the molecular recoil frequency due to the formation of bound pair state, which dominates the response 
tie response function S^, given in Eqs. fl2.34 2.43L receives contribution from non-interaction terms 
(2.351, which peaks at the atomic recoil frequency. On the BCS side, the interaction is too weak so the response 
is dominated by the non-interacting terms. On the BEC side, the peak of response function moves toward the 
molecular recoil frequency as Sf±. 



T/T F =2 



T/T F =3 



0.2 - 



I I I I I I I I I I I I I I I I 



3, 

<T? 0.2 



S^(q,co) 



(b)l/k F a=0 



CZ3 



S ti (q,co) 

I I I I I I I I I I I I I I I I I I I 



0.2 



(c)l/k a=-1.0 



J I I I I I I I I I I I I I I I I I L 



~i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 



(d)l/k a=1.0 




- 0.2 



- 



S^(q,co) 



(e)l/k F a=0 



- 0.2 



S n (q,(0) 



I I I I | I I I I | I I I I | I I I I 

(f)l/k„a=-1.0 



0.2 



j i i i i i i i i I i i i i i i i i i_ 



1 

00/ 00 



R 



1 

oo/oo 



R 



FIG. 1: (Color online) Sff(q = 3&;f,^) and Sfi(q = 3kF,oj) as function of frequency at different 1/kpa and at T — 22> (left) 
and T = 37> (right). 



B. Interaction dependence of dynamic response function 



In Figure [2] the dynamic density (left) and spin (right) responses at T = 3Tp, obtained by Eqs. (2.1 2.2), are 



shown for various interaction strengths. The density response function peaks at atomic recoil frequency on BCS side. 
As interaction increases, the peak becomes red-shifted. On BEC side, the density response becomes peaked around 
molecular frequency. These features qualitatively agree with recent Bragg spectroscopy experiment [2J, although the 
latter experiment was performed at extremely low temperature in the superfluid phase. The spin response has broad 
peak around atomic frequency both on BEC and BCS sides, which again agrees with the experiment [2J. 




FIG. 2: (Color online) Interaction dependence of dynamic density (left) and spin (right) responses at T = 3Tf. 



C. Temperature dependence of dynamic response function 

In Figure |3j the dynamic density and spin responses at unitarity are shown at different temperatures. The density 
response has a relatively narrow peak around ujr/2 at T = 2Tp, and becomes broadened and shifted to higher energy 
as temperature increases when non-interacting term becomes increasingly more important. On the other hand, the 
spin response peaks around ojr. 



D. Static response function 



The static response functions are obtained via following integrals, 

/oo 
dioS a<J >{q,Lo) 
-OO 



(3.1) 



In FigurcHthc static response functions Sff (q=3kp) and Sf± (q=3fcp) are shown as function of 1/kpa at T = 3Tp. 
Sff (q=3kp) increases from about one in deep BEC side when 1/kpa reduces, and it decreases quickly approaching 
unitarity. The rise and fall is due to contribution from transition between bound state and scattering state on the 
BEC side in Eq. (2.36), which vanishes at unitarity. On the other hand, S-j-j, (q=3kp) changes more smoothly from 
BEC to BCS side. Note at unitarity when T = 3Tp, the fugacity is about 0.13, therefore one expects the higher order 
terms contribute at about 10-20% of second order term displayed here. 
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FIG. 3: (Color online) Temperature dependence of dynamic density and spin responses at unitarity. 

IV. CONCLUSIONS 

In this work, the dynamic response functions of strongly interacting Fermi gas in homogeneous space are investigated 
in a virial expansion to second order. The dynamic density response function is found to exhibit transition from atomic 
response to molecular response, as the interaction strength increases and the system undergoes BCS-BEC crossover. 
The spin response function has a broad peak around atomic recoil frequency at high temperature. These features 
qualitatively agree with recent experiments via Bragg spectroscopy. The virial response is exact at low density and 
high temperature, when fugacity is a small number. The fugacity discussed in this work is relatively small z < 0.25. 
Therefore it provides a benchmark for the many-body response functions. 

Qualitatively, the response functions of strongly interacting fermi gas in homegeneous space show similar charac- 
teristics as those in trapped fermi gas. In fact they may be related by a local density approximation, which would be 
interesting to study in future work. 

To second order virial expansion, the dynamic response function in homogeneous space can be explicitly written 
down in compact and closed form with simple integrals. This illustrates that the virial expansion to response function 
may be also applied to other strongly interacting many-body system at low density and high temperature, like neutron 
matter and nuclear matter in supernova. 

It is a pleasure to acknowledge the encouragement from Joe Carlson, Chuck Horowitz, and Sanjay Reddy, as well as 
discussions with C.-C. Chien, J. Drut, S. Gandolfi, and Y. Nishida. The work is supported by a grant from the DOE 
under contract DE-AC52-06NA25396 and the DOE topical collaboration to study "Neutrinos and nucleosynthesis in 
hot and dense matter" . 
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FIG. 4: (Color online) Static response functions Sff (q=3kp) and Sfi (q=3kp) as function of 1/fcpa at T = 3Tp. 

Appendix A: S-wave pair states in a <5 3 (r) potential 
The S wave scattering state is known to be, 



2 sin(pr) -pacos(pr) / 2 sin(pr + rj p ) 

<Ps(pr) = \ 7z 2 2 , >oo = \ >oo, (Al) 

7r(l + p z a z ) pr \ tt pr 



when rj p equals zero, it reduces to non-interacting S wave state. To be precise, there should be a factor — sign(a). 
However this overall phase does not influence the results in the paper. The bound state on the BEC side (a > 0) is, 

1 e~ r l a 

Mr) = -j= , (A2) 

with the binding energy, 

E b = V (A3) 
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Appendix B: Sum of product of two Legendre Polynomials 

The infinite sum over product of two Legendre polynomials can be carried out utilizing following relations, 

^(-1)'(2Z + 1)P^) 2 = ^), XG[-1,1], (Bl) 
z>o 

£)(2Z + l)fl(aOQ,(aO = 0, (B2) 
;>o 

X)(-l)'(2i + l)fl(x)Q,(a:) = ix e [-1, 1], (B3) 

2>0 

£(2Z + l)Q,(x) 2 = (|) 2 £(2Z + l)fl(*) 2 + -^— r (B4) 

l>0 l>0 

£(-l)<(2Z + ljQ^) 2 = ~ log |i±||. (B5) 
i>o x 



Note in Eq. (B4) the two sums are separately divergent when x € (—1, 1), but their difference is finite. 



Appendix C: Virial expansion of density 

On BCS side, the following equation for number density, 

" x ' - <"^> s,! - 2 r it * + ^(WMi-li, (CD 



2 67r 2 V^Wo i + e*/z \ fk F a 
is used to solve fugacity z. 

On BEC side, the equation for number density is changed accordingly due to the additional contribution from 
bound state, 

"4 = {i ''T -±r*rAr + ^< f '*">' ( > - a**/! ^) < C2 > 

2 6?r 2 v/ttJo l + e t /z \ l \ T k F a' y ' 
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